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An experimental research concerning highly underexpanded jets made of different gases from the 
surrounding ambient is here described. By selecting different species of gases, it was possible to vary 
the jet-to-ambient density ratio in the 0.04 to 12 range and observe its effect on the jet morphology. 
By adjusting the stagnation and ambient pressures, it has been possible to select the Mach number 
of the jets, independently from the density ratio. Each jet is therefore characterized by its maximum 
Mach number, ranging from 10 to 50. The Reynolds number range of the nozzle is 10^ to 5 • 10**. 
The spatial evolution of the jets was observed over a much larger scale than the nozzle diameter. 
The gas densities were evaluated from the light emission induced by an electron beam and the gas 
concentrations were obtained by analyzing the color of the emitted light. The results have shown 
that the morphology of the jets depends to a greater extent on the density ratio. Jets that are 
lighter than the ambient exhibit a more intense jet-ambient mixing than jets that are heavier than 
the ambient, while the effects of changing the jet Mach number do not seem to be too large in the 
explored range. These results can be expressed by means of two simple scaling laws relevant to the 
near field (pre-Mach-disk) and the mid-long term field (post-Mach-disk), respectively. 

Keywords: underexpanded jets, mixing layer growth, density and concentration gradient, barrel shock, Mach 
disk 



I. INTRODUCTION 

The experiment here described has been devoted to the 
study of hypersonic underexpanded jets. In the last few 
decades, these jets have received a great deal of attention 
because of their importance in several fields such as basic 
fluid dynamics and astrophysics, and in applications for 
the aeronautical and mechanical industries. 

The near field of underexpanded supersonic jets is 
characterized by a well known structure constituted by 
an intense expansion up to a high Mach number which 
ends in a normal shock, the Mach-disk, and which is sur- 
rounded by a barrel-shock. This structure was described 
in the early theoretical and experimental works by Ashke- 
nas and Sherman [l|, Crist et al. ^. Later on, the high 
expansion present in these jets was represented by mod- 
els that linked the Mach numbers and shock properties 
to the entropy change of the fluid from stagnation to the 
final after-disk state Q . 

Very high pressure ratios can also give rise, at the 
end of the expansion, to a transition from a continuum 
to an almost coUisionless flow. The relevant nonequi- 
librium phenomena have been the subject of numerous 
investigations, see the work of Cattolica et al. [J|, who 
made kinetic temperature measurements and qauntified 
the anisotropic partition of the random thermal energy 
along the direction parallel and normal to the mean flow. 

The properties of underexpanded jets have also been 
considered in astrophysics, where the evolution of jets 
observed in young stellar objects is an important issue 



[5| . In this field, some information may emerge from the 
study of the underexpanded jet morphology as a function 
of the Mach number and jet-to-ambient density ratio |6|. 

Numerous investigations have focused on underex- 
panded jets, and there is an extensive literature on the 
subject which includes compared numerical simulations 
and laboratory experiments (for example [7|, |8|), both in 
basic and applied research. In general, the near jet is 
described in a detailed way, and complete information is 
given on its velocity, pressure, temperature and density 
field [3, [l3l : but most of the literature on these fiows does 
not report data about the subsequent part of the jet, al- 
though there are a few exceptions [5|, Q. Moreover, it 
is not easy to find data about the dependence of the jet 
properties on the jet-to-ambient density ratio. 

The experimental facility employed here is suitable for 
investigating the spatial evolution of a jet over a scale 
of nearly two orders of magnitude larger than the for- 
mation scale, which is here defined as the nozzle throat 
diameter. This makes it possible to investigate a long 
region downstream from the Mach disk. Furthermore, a 
unique feature of the present experiment is the possibil- 
ity of studying the effects of two flow control parameters, 
the stagnation-to-ambicnt pressure ratio and the jet-to- 
ambient density ratio, independently of each other. The 
measurement techniques are based on the electron beam 
technique. The fluorescence induced in plane sections of 
the flow field are acquired as digital images and processed 
by means of a specially developed algorithm, which yields 
density and species concentration maps |ll| . 



A number of jets has been studied here by varying the 
stagnation-to-ambient pressure ratio (which is the coni- 
pressibihty parameter that sets the maximum Mach num- 
ber and the Mach disk location in the near jet) and the 
jet-to-ambient density ratio over a range. The variation 
of this last paramenter was possible due to the choice of 
different pairs of gases from a set of a few noble gases, 
such as Helium, Argon or Xenon (this is the parameter 
that sets the transverse density gradient in the interme- 
diate and far jet). The density and concentration fields 
were determined for each jet, and the relevant Reynolds 
number was carefully characterized. It has been shown 
that the ambient density conditions play an important 
role on jet dynamics and morphology. The obtained 
results could help us to understand their stability and 
transport properties. In particular, it has been shown 
that light or underdense jets (where the jet density is less 
than that of the ambient) spread out faster than heavy 
or overdense jets. This is in agreement with previous 
findings relevant to finite-thickness one stream compress- 
ible shear layers with uniform density (Gropengiesser, 
1970) [i2l and to compressible turbulent shear layers (Pa- 
pamoschou and Roshko, 1988) [IJl. In the present study, 
we have been able to link the spreading angle before (near 
field) and after (mid-long term field) the Mach-disk and 
the jet/ambient density ratio by means of two simple 
scaling laws. 

This work is organized as follows: section II presents 
the flows under study and the laboratory facilities, sec- 
tion III describes the data analysis methods, section IV 
contains the experimental results and section V the rele- 
vant discussion. 



II. JETS UNDER STUDY AND LABORATORY 
SET-UP 



CCD camera 




FIG. 1: Experimental setup. The ambient gas inlet is 
mounted onto the vacuum vessel wall. The nozzle section 
is shown in the lower right corner. 



TABLE L Capabilities of the experimental system. 



Vessel length 




2m to 5m (modular) 


Vessel diameter 




0.5m 


Vessel (ambient) pressure 


Pamh 


1.5 to 200 Pa @ 300 K 


Nozzle throat diameter 


D 


2mm 


Jet stagnation pressure 
Stagnation/ambient 
pressure ratio 


PO 
Po/Pamb 


2000 to 2 ■ lO^Pa 
10 to 1.3 ■ 10^ 


Electron gun voltage 




up to 20kV (16kV typ.) 


Electron gun current 




up to 2mA (1.5mA typ.) 



The present experiments are carried out in a system 
that was specifically designed for the study of hypersonic 
jets [a, [Tl|. The system facilities consist of a cylindri- 
cal vacuum vessel equipped with suitable nozzles, and an 
electron gun that makes it possible to obtain the fluores- 
cent plane sections of the jets, which have here acquired 
as digital images using a CCD color camera. The setup 
inside the vessel is sketched in fig. [1] 

Here the tested fiows are strongly underexpanded jets 
obtained from a truncated sonic nozzle that flow along 
the longitudinal axis of the vacuum vessel. The vessel has 
a diameter of 0.5 m which is always 5 times larger than 
the diameter of the jets, which ranges from 0.01 to 0.1 m. 
The vessel is equipped with 2 valve systems that control 
the jets issued by the nozzle and the ambient gas; this 
offers the opportunity of using different gases for the jets 
and the surrounding ambient. The jet stagnation pres- 
sure pq and the ambient pressure Pamb are determined by 
adjusting the opening times of the 2 valve systems, using 
by a suitable electronic device. The essential data of the 
experimental system are shown in table HI 



The accuracy of pressure transducers is 1% on the 
Pq measurements and 0.25% on the Pamb measurements. 
The upper ambient pressure limit is essentially set by ac- 
cording to the proper conditions for electron gun work- 
ing, i.e. pressure below 200 ^ 300Pa and density number 
no greater than lO'^'^m"^ in order to avoid excessive ab- 
sorption of the electrons by the gases and light intensity 
saturation. The camera has a high sensitivity, IMpixel 
CCD, with a Bayer RGB filter, the operating mode is 12 
frames/s for a 83 ms exposure or 24 frames/s for a 41 ms 
exposure. 

The flow time scale, defined as the ratio between jet 
diameter and local speed of sound, ranges from 10/iS in 
the nozzle throat to 1ms at the end of the jet expansion, 
upstream from the Mach disk. The typical jet outflow 
time is 0.5s, that is, hundreds or thousands of time scales, 
and changes in the ambient can be observed in this lapse 
of time due to the jet/ambient mixing. However, due 
to the large mass flow of the primary pumps, fast varia- 
tions in the ambient during the outflow time are avoided, 
and since the image exposure is sufficiently shorter than 



the outflow time, the large-scale properties of the jets 
show only very small variations over a single exposure, 
therefore the jets can be considered quasi-steady in each 
image. 

Choosing the jet/ambient pairs from Helium, Argon 
and Xenon and varying pq and Pa in the allowed ranges 
permits variations of the jct-to-ambicnt density ratio 77 = 
Pjet/ Pamb, ranging from 0.04 to 12 in the near jet zone 
and from 0.2 to 4 in the intermediate and far jet zone, to 
be obtained. 

In these experiments, the Mach number of the jets in 
the near zone is very high, since the jet gas may ap- 
proach the limit velocity just upstream from the Mach 
disk, therefore M values of up to 50 are possible. The 
Knudsen number K is always in the continuum regime, 
except for the particular case of helium jets with very 
high po/pamb ratios, where K may approach 0.5 in a 
small region upstream from the Mach disk. It should 
be noticed that in this case and in general at the end 
of the expansion, the accuracy of the measurements can 
deteriorate [ll[. In the following sections, the reported 



accuracies always take this problem into account. 



III. DATA ANALYSIS 

A. Density and concentration measurements 

Here the measurement method described by Belan et 
al. |ll| is used to obtain density and concentration maps 
of jet and ambient gases from the original digital color 
images. The method is based on the analysis of the flu- 
orescent emission / from a gas excited by an electron 
beam, which at the low pressures considered here can be 
expressed in the form 



(1) 



where n is the gas number density and k a proportionality 
coefficient. 

The limits of validity of this approximated law are 
stated by the following conditions: 

- low density (n < lO^^m"'^), which ensures linearity and, 
for gas mixtures, a decoupled emission (i.e. the total 
emission is simply the sum of the partial intensities due 
to different species) 

- bounded range for the gas temperature {70K < T < 
900K), which cannot be so high as to have a visible emis- 
sion, or so low as to have an emission affected by nonequi- 
librium phenomena. 

These conditions have been widely satisfied by the 
tested jets, and in particular in the jet/ambient mixing 
zone. The only exceptions are: 

- a small region upstream from the Mach disk, where the 
temperature may fall below 70K and the uncertainties on 
the measurements can grow by more than 15%. 

- a small region close to the nozzle, where the density is 
too high and the emission is saturated. 



The operations performed by the image post- 
processing algorithm are summarized hereafter, while the 
relevant details are reported in appendix [A) The jet gas 
concentration along an image, Zj^t, is obtained reading 
the RGB values of each pixel that form the color ratio 



Ci aiR + biG + ciB 



C2 



a2R + b2G + C2B 



(2) 



where the six coefficients ai,bi, Ci are chosen on the basis 
of the gas pair involved in the experiment, and using the 
relation 



Zjet 



^amb\ 



■ib^r 



'^arnbi '^jeti ~r y'^jet2 '^amb2 ) ^ 



(3) 



where the four coefficients kambi,kamb2,kjeti,kjet2 are 
known from calibrations performed on pure gases. The 
concentrations of the two gases are obviously linked by 
Zamb = 1 — Zjet ■ Oucc the Concentrations are known, the 
density is given by the equation 



P = Ci 






(4) 



where mamb,iTT'jet are the molecular masses of the gases 
in the ambient and in the jet, respectively. 



B. Measurements on plumes and mixing layers 

The aim of this section is to define suitable criteria 
to compare the different jets under study and to exploit 
their dependence on the control parameters. To this aim, 
the measurement method explained above will be used 
to determine the jet properties by means of density and 
concentration values. 

At this point, it is important to consider the under- 
expanded jet structure in detail, referring to fig. [2j In 
the short range upstream from the Mach disk, which is 
located at a distance L from the nozzle (near jet or barrel 
zone), the Mach number increases to its maximum value 
Mmax just upstream from the normal shock, whereas the 
density ratio decreases very rapidly: this happens both 
for light and heavy jets. In the present case, the jet prop- 
erties essentially depend on the stagnation-to-ambient 
pressure ratio. 

The typical trends of radial the concentration and den- 
sity curves in this zone are shown in figure [3^. It should 
be noticed that Zjet does not vanish in the outer zone: in 
fact, as the testing time increases in these experiments, 
the mass fiow of the jet gas, which is not completely 
balanced by the primary pumps, despite their power, is 
partially mixed with the environment. In a given ex- 
periment, for example Ar in He, the Ar jet is therefore 
found to travel in an environment of Ar mixed with He. 
Thus, in these experiments, the ambient contains an in- 
eliminable part of the jet gas, but this is not a limit to 
the present research since this effect is directly taken into 
account in the density and concentration measurements. 




Near Jet 



Intermediate Jet 



Far Jet 



FIG. 2: Structure of a strongly underexpanded jet. L is the nozzle-Mach disk distance or barrel length, a is the plume 
growth angle measured on an arc of length I, r is the conventional jet radius, and /3 is the mixing layer spreading angle in the 
intermediate jet. 
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FIG. 3: Jet cross-sections in the near (a) and intermediate (b) regions. 



Here, in the near jet, a 'mixing layer width' or 'plume 
width' can be defined as the radial distance between the 
barrel shock and the outer jet boundary (figure [3]): the 
radial position of the barrel shock, yin , can be identified 
as the point where the density curve has the maximum 
slope, 



d"p 
d^ 







(5) 



and the outer jet boundary can be identified by condition 



■^jet 



dy' 







(6) 



Vout 



i.e. where the derivative dzjet/dy, with respect to the ra- 
dial coordinate y, reaches its maximum value. It is easy 
to see that this is analogous to the classic 'half width at 



half height' criterion that is extensively used for jets in lit- 
The y value where condition 



crature 



14 



is satisfied 



can also be taken as a conventional jet radius r = yout, 
bearing in mind that the jet width 2youti obtained as 
the distance between opposite boundaries, has a relative 
and not absolute meaning (for example, other quantities 
proportional to yout could be taken as conventional jet 
radii without any loss of generality when similar jets are 
compared). 

Therefore, in a given radial section, the plume width 
5mix can be defined as 



^mix — Vout yin- 



(7) 



For the sake of comparison, the jet cross section that 
should be considered when measuring Smix should be con- 
ventionally taken at the same axial position X in each ex- 
periment, here, X = 0.8L, referring to the nozzle-Mach 



disk distance L (barrel length), has been chosen; X is 
not too close to the Mach disk zone to make the mea- 
surements too difficult. Once 5mix has been determined, 
it is possible to calculate the plume growth angle a: 



Reynolds number is used for this purpose, and is defined 
in the barrel zone {x < L) as: 



a ^ eiTctan{5mix/l) 



(8) 



where I is the relevant arc length, that is the length of 
the barrel shock between the nozzle and the measurement 
cross section for Smix- 

In the long range (intermediate and far jet), down- 
stream from the Mach disk, the Mach number remains 
quite low, and the jet properties essentially depend on 
the density ratio. There are also density fluctuations in 
some underdense jets, which are shown in figure [2] as thin 
line beams. Here, the mixing layer grows both inwards 
and outwards: it is possible to use the same definition 
as above for the outer boundary of jet and mixing layer, 
while a criterion based on concentration and density val- 
ues, even though technically possible, results to be too 
noise-sensitive for the inner boundary. Nevertheless, the 
long scale behaviour of the jets can be described by a 
spreading angle /3, defined in terms of the outer bound- 
ary position. This angle represents the increase in the 
mixing region in the intermediate field where an annular 
mixing layer surrounds the jet core; in the final far field, 
where the jet is fully developed and the flow becomes un- 
compressible, /3 instead represents the spreading angle of 
the whole jet, but this field is beyond the scope of this 
work and will not be considered. 

In order to identify the jet boundary in the intermedi- 
ate field only the cross-sectional Zj^t curve is required; a 
typical trend is shown in figure |3]3. After choosing a pair 
of axial positions xi,a;2 in the range x > L, the outer 
boundary is found using criterion ([5]), and two values of 
the jet radius ri = yout{xi) and r2 = yout{x2) are ob- 
tained. Then, the spreading angle is obtained from the 
increments Ar = T2 — ri and Aa; = X2 — xi through 
relation 



/3 = arctan(Ar/Aa;) 



(9) 



These measurements are usually performed choosing sev- 
eral xi, xi pairs and calculating /3 as a mean value. The 
X\^X2 pairs are choosen in the x > 1.2L range to avoid 
setup effects on the mixing layer. 

In order to ensure the proper evaluation of the concen- 
tration and density curves, together with their deriva- 
tives, all the images analyzed in this work have been fil- 
tered to improve the signal-to-noisc ratio by means of 
standard local averages in the small neighborhoods of 
each pixel. The non-uniformity of the electron sheet in 
the images is also compensated for. 



C. Reynolds number 

The description of the tested jets should be completed 
with their Reynolds numbers. Here, the longitudinal 



RBx 



Uxp 



(10) 



where U ~ Ale — M (x) y^ jTZT (x) is the ccnterline ve- 
locity, the local Mach number M{x) is calculated accord- 
ing to Ashkenas et al. and Young [l|, Q and T{x) = 
T[M{x),Tq] is the isoentropic temperature relevant to 
the stagnation temperature Tq and to M(x). The den- 
sity p is measured directly on the images, and fi{T) is 
calculated for each gas by means of the Sutherland law 
where possible and of suitable models at the lowest tem- 
peratures J16l - [l9j . For the sake of comparison among 
different jets, in what follows for each jet it will be re- 
ported the Reynolds number attained just upstream of 
the normal shock (x ~ L). 

A remarkable property of these flows is the range of 
Reynolds numbers covered by each jet; actually in these 
jets the Reynolds number, whether diameter based or jet- 
length based, falls off suddenly after the Mach disk. This 
variation can be of several orders of magnitude, since 
the strong temperature rise increases the viscosity dizzy, 
while the velocity decrease is smaller and the space scale 
varies slowly. The same problem arises if the Reynolds 
number is estimated for the mixing layer zone. An exam- 
ple can be given for a typical Argon jet having a pressure 
ratio Po/Pamb = 10*: in this case, the diameter based 
Reynolds number is 2.4 ■ 10* in the throat, then it in- 
creases up to 6.7 • 10^ upstream of the Mach disk, and 
falls off to about 83 downstream. 



IV. RESULTS 

The present results have been obtained considering 280 
image sequences, each consisting of 18 or 36 images of 
jets under quasi-steady conditions; a wide range of Mach 
numbers and density ratios, ranging from underdense jets 
to overdense jets, is therefore available. Among these 
images, those with the best S/N ratios were selected to 
obtain the presented results. In what follows, the im- 
age processing techniques are first applied to a selected 
sample image as an example; the general results are then 
presented. 



A. An image processing example 

Figure 2] shows an Xenon jet flowing in a Helium ambi- 
ent. The jet stagnation pressure is pq ~ (4.10±0.20) • 10* 
Pa, while the pressure of the Helium ambient is Pamb = 
(9.86 ±0.1) Pa. The flgure is overcontrasted for the sake 
of clarity. The measurements have been performed along 
the axis and on different cross-sections (see below). 

Figures [5] and [6] show examples of the axial and cross- 
sectional density and concentration measurements per- 
formed in flg. m Each figure contains sample error bars. 




100Q- 



FIG. 4: Example: Xenon jet in a Helium medium. Pressure 



ratio pa/pamh ~ 4.3- 10 , Mach before the normal shock • 
Reynolds number at the nozzle exit = 2.5 • 10*. 



'41, 



as typical values of the uncertainties, which are obtained 
by propagation of the errors, accounting for noise values 
along the image and for calibration coefficient uncertain- 
ties in the image processing algorithms. It can be seen, 
in fig. [5] that the density values along the jet axis agree 
well with the relevant theoretical isoentropic expansion, 
a fact which should presumably be related to the heavy 
weight of the Xenon atoms. However, this is not a gen- 
eral property of strongly underexpanded jets, since the 
gap between the measured and isoentropically calculated 
values for undcrdcnse jets, for example He in Ar, is larger 
[Il|. The Mach disk position is in good agreement with 
Young's theoretical estimation [31 . 

Three radial density curves can be seen in fig. [5] at 
different axial stations, namely x — 32D ~ 0.8L, x = 
56D = lAL and x = 801? = 2L, in terms of the initial 
diameter and the barrel length. Here the uncertainties 
are in general of the order of ±10% of the local axial 
density value, except in the first x/ D = 32 curve, where 
the axial zone has a lower S/N ratio, which is usual in a 
dark underexpanded zone. The barrel shock position is 
clearly visible in the curve for x/D = 32, the other two 
curves show the effects of species mixing in the considered 
range. 

Figure [7] shows the Xe concentration map relevant to 
figure 131 Here, the almost pure Xe zone [zjet > 0.9) 
slowly becomes narrower as the jet travels along the ax- 
ial direction. There is a visible spreading in the other 
contour lines, but the width of the region between the 
maximum variations in the Zjet curve does not change 
fast, which means the jet boundaries spread out quite 
slowly according to rule ^. However, in the following 
results it will be seen that the spreading of underdense 
jets is instead remarkably larger. 
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FIG. 5: Axial density for the Xenon jet visualized in fig. |4l 
A sample error bar is shown on the right. The zone on the 
left cannot be analyzed because of the high density values 
(saturated output). The dashed line shows, as a reference, an 
equivalent isentropic expansion followed by a normal shock. 
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FIG. 6: Cross-sectional density curves for the Xenon jet vi- 
sualized in fig. [4] Sample error bars are shown on both sides 
of each curve. 




10 20 30 40 50 60 70 80 90 
x/D 



FIG. 7: Concentration map of the Xenon jet visualized in fig 
m The zone marked with "N" cannot be analyzed, since the 
relevant pixels are in the density saturation range, fall outside 
the electron sheet, or are influenced by the nozzle image. 



B. Plume spreading in the near jet 

This section deals with results referring to the region 
between the nozzle and the Mach disk. Among the 
numerous available jets, the ones selected according to 
the best S/N ratio are: 

8 Helium jets traveling in the Xenon ambient (labeled 
He/Xe) 

13 Helium jets traveling in the Argon ambient (labeled 
He/Ar) 

14 Argon jets traveling in the Helium ambient (labeled 
Ar/He) 

8 Xenon jets traveling in the Helium ambient (labeled Xe/He) 

Since the jet stream is accelerated in this zone, each 
jet will be characterized by the Mach number Mmax at 
the end of the expansion, upstream from the Mach disk, 
for comparison purposes among the jets. M^ax is calcu- 
lated according to Young 0]. In an analogous way, the 
maximum longitudinal Reynolds number Rcx, calculated 
as in ijIIICI bv equation ([T0|) at the end of the expansion, 
wiU be reported for each jet. 
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FIG. 8: Parameter domain Mmax rj (maximum Mach number 
- density ratio) of the tested jets. Sample error bars are shown 
on the right. 

Since the independent parameters in these experiments 
are the density ratio and the Mach number, each se- 
lected jet corresponds to a point in the parameter do- 
main Mmax "Hi as shown in figure [5] Some sample values 
of the uncertainties are shown on the right side of each 
curve, since the parameter pairs Mmax V are not set di- 
rectly because of to the experimental setup; their values 
depend on the pressure ratios po/Pamb and the gas pairs, 
which are instead set directly and which permit the un- 
certainties of Mmax and 77 to be calculated. The reported 
density ratios 77 refer to a cross-section lying at 0.8L from 
the nozzle, that is, a value which is close to but not inter- 
cepting the Mach disk. The density ratio, on each cross- 
section, is obtained by measuring the maximum pjet, i-e. 
the jet density at the outer surface of the barrel shock 
(the p-curve peak in fig. |3^) and measuring pamb in the 



outer far field. It is also possible to determine ry taking 
Pjet on the jet axis, but this is less diagnostic, because 
the strong expansion in this zone makes the density ratio 
less dependent on the gas species, and makes the S/N 
ratio worse. For example, the overdense jet in figure [6] 

gives Pjet, max/ Pamb = 4.83, whilc Pjet.axis / Pa,nb = 1-35 

at the same cross-section. 

All the jets considered here have nozzle Reynolds num- 
bers, based on the throat diameter, ranging from 10^ to 
5-10^. For the selected jets, the maximum Reynolds num- 
bers attained upstream from the Mach disk are reported 
in fig. ini as functions of the maximum Mach numbers. 
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FIG. 9: Maximum longitudinal Reynolds number Re^ up- 
stream from the Mach disk as a function of Mmax- Sample 
error bars are shown on the right 

Figure [10] reports the growth angles of the jet plumes 
as functions of Mmax for different gas pairs, ranging over 
several density ratios. The relevant r\ values, reported 
in the caption, are the same as those in figure [51 The 
plume growth angles a arc obtained from formula ([8|) , as 
explained in §111 Bl The ucertainties are calculated by 
Gaussian propagation as above. 

The general trend exhibited by these results is that 
underdense jets spread out more than overdense jets do. 
The plume spreading angle a seems to depend weakly 
on the Mach number Mmax- it decreases slowly for un- 
derdense jets as Mmax grows, whilst it becomes nearly 
independent for overdense jets. The Reynolds depen- 
dence of a is also analogous, i.e. a diminishes slowly for 
underdense jets as Rex grows. Apart from these prop- 
erties, it is worth noting that the density gradients in 
overdense jets are sharper , and it can be shown that the 
shocks are also thinner, a result that has been confirmed 
by acquiring monochromatic intensified images with very 
short exposures [l5|. 

The global trend of the a values can be expressed by 
means of a scaling law in the form 



« = a + 6A/("+'='''; 



(11) 



this law best fits the experimental data with the coef- 
ficients a = 4.67, h ^ 36.9, n ^ -0.235, k = -0.315, 
which give an estimated error variance of 0.627. The 
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FIG. 10: Spreading angles of the plumes outside the barrel 
shock vs Mmax for several jet/ambient gas pairs The density 
ratio ranges are 0.50 < -q < 0.62 for He/Xe, 0.95 < t? < 1.20 
for He/Ar, 1.54 < r) < 2.40 for Ar/He and 3.40 < 77 < 12.0 
for Xe/He. 




FIG. 12: Spreading angles of intermediate jets vs density ratio 
(the axial Mach number downstream from the Mach disk is 
about 0.45 for all the jets). The green line is a fit according 
to Papamoschou and Roshko (1988) (the red dashed line is 
the best scaling law fit, see formula (|13p. 



relevant surface a{Mjnax , v) > superimposed onto the ex- 
perimental data, is shown in figure 1111 
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FIG. 11: Scaling law (|ll|l for the plume spreading angle as a 
function of the Mach number Mmax and density ratio rj. 



C. Mixing layer spreading in the intermediate jet 

This section contains results referring to the interme- 
diate jet region, downstream from the Mach disk. The 



mixing layer spreading angles /3 for the intermediate jets 
are obtained from formula ([9]), as explained in ijIIIBI 
These angles are determined from data collected on an 
axial domain 1.2L < x < Xmax, where Xmax ranges from 
2.5L to 5i, depending on the jet under test, and avoid- 
ing positions where the mixing region reaches the jet axis. 
Figure [T^] reports the /3 values as a function of the den- 
sity ratio, where the very weak dependence on the Mach 
number has been omitted (see below). The r] values are 
determined from the axial mean value of the jet center- 
line density, which varies slowly in the considered range, 
and from the density value of the surrounding far ambi- 
ent. These measurements can be very sensitive to image 
noise, depending on the involved gases, and require a 
wide image field: for these reasons, figure [T^] just reports 
a subset of the available jets. The uncertainties have been 
calculated with the same method as above; the resulting 
values arc not small. 

Here the dependency on the Mach number is weak: 
the Mach number Mi on the jet axis after the normal 
shock is in fact always subsonic, and very close to the 
minimum theoretical value (.448 < Mi < .455), since 
the flow upstram of the normal shock has Mmax > 10 
for all the jets studied here. After the normal shock, the 
jet pressure is also very close to the ambient pressure. 
Thus, the compressibility effects in this region seem to 
be small and the jet evolution is essentially determined 
by the density ratio. 

Nevertheless, although the jet core becomes subsonic 
and nearly incompressible after the Mach disk, an annu- 
lar supersonic region in the outer jet may survive over a 
short distance, which is less than 2L in these experiments. 
This is due to the different behaviour of the inner and 
outer streamlines: the inner streamlines fall from Mmax 
to Ml ~ 0.45 after crossing the Mach disk, whereas the 
outer streamlines that cross the barrel and the reflected 
shock visible in fig. [2] - both obliques - have a larger 



exit Mach number than unity (see, for instance, Norman 
and Winkler [5|). Here the flow remains supersonic until 
the lateral mixing region expands over the entire jet sec- 
tion, and this can give rise to further expansions along 
the je t up to the sonic condition, followed by recompres- 
sions |20| . These phenomena were in fact observed in the 
undcrdcnsc Hc/Ar and Hc/Xc jets at pressure ratios of 

100 < po/parab < 1000. 

On the other hand, if the effects of the supersonic 
tongues and the possible velocity oscillations on the axis, 
in the intermediate jet zone are disregarded, the flow can 
be considered as a core stream surrounded by an annular 
mixing layer, with the static pressures of the core and 
ambient very close to each other. Under these hypothe- 
ses, and expressing the mixing layer thickness measured 
along the x axis with respect to a given origin xq, as 5{x), 
the growth rate 



5'{x) 



X — Xq 



can be estimated by the measured angle /3, since 2/3 ■^ 
arctaufJ' ^ 6', under the hypothesis of small (3. The 
measured angles can be compared with the Papamoschou 
and Roshko model [13[, which can here be expressed by 
a law of the kind 



/3 = c[l + r;-i/2], 



(12) 



where c is a constant. The best flt based on this law 
has a coefficient c = 3.01, see formulas 14 and 15 in [13l |. 
The resulting curve is shown in figure [12] together with 
the present experimental data. 

The general trend in the intermediate and far regions 
may now be outlined as follows. In the undcrdcnsc jets, 
the mixing layer spreads fast, even toward the jet axis. 
Visualization and measurement attempts have shown 
that the mixing layer growth in these jets leads to a fast 
destruction of the jet core, typically on a length of the 
order of ^ AL. These jets also have weaker gradients. 
Instead, in the overdense jets, the mixing layer spreads 
slowly, therefore the jet core has a longer path, usually 
> 4i. 

The global trend of the /3 values can also be expressed 
by means of a general scaling law in the form 



(3 = a + brf' 



(13) 



which has the best fit on the experimental data with the 
coefficients a = 1.54, h = 4.61, n = —0.32, while the 
estimated error variance is 0.482. The relevant curve 
/3(?7), superimposed onto the experimental data, is shown 
in figure [T^l and results to be very close to model (IT^ , 
except for very small values of r/. 



V. CONCLUSIONS 

The present experiment has mainly focused on the ef- 
fects of the jet-to-ambicnt density ratio on the morphol- 
ogy of strongly undcrexpanded jets, on both short and 



long axial scales, in the range where compressibility is 
important. All the tested gases are monoatomic, and 
therefore there is no dependence of the stream proper- 
ties on the specific heats ratio 7, which has simplified 
the characterization of the density effects. The Reynolds 
number in the nozzle, whihc is diameter based, ranges 
from lO'^ to 5 ■ 10"*, whereas the Reynolds number in the 
expansion, based on the axial length, ranges from 2-10'^ 
to 5 • 10* for light jets and from -10^ to 2 ■ 10^ for heavy 
jets in the same Mach number interval, due to the differ- 
ent densities involved. All the spreading angle measure- 
ments have been based on gas density and concentration 
values; the resulting mixing layer spreading dependence 
on the density ratios show that the jet/ambient mixing 
is in general more efficient in underdense jets. In partic- 
ular, it has been shown that the plumes in the near field 
spread up to three times wider in underdense jets than in 
overdense ones; in the intermediate field considered here, 
which is subsonic downstream from the Mach disk, the 
spreading of light jets is up to two times wider than for 
heavy jets. This different behaviour is slightly influenced 
by the Mach number in the near field, where it has been 
shown that increasing Mach numbers have a weak colli- 
mating effect, and reduce the plume width. The whole 
body of results has been synthesized by means of two 
simple scaling laws, which are valid in the near and in 
the intermediate fields, respectively. 



Appendix A: 



This section reports details of the image post 



processing algorithm used in this work, according to [11 



The method is based on equation ([T]), i.e. the linear law 
I = kn for the fluorescent emission / of the excited gas 
where n is the gas number density and k a proportional- 
ity coefficient. 

The emission / is in fact a spectral superposition, 
therefore equation ([T|) also holds for the three colors (R, 
G, B) acquired by each pixel of the camera sensor and 
for any of their linear superposition: 



R = 


= kRU, 




(Al) 


G . 


= ken, 




(A2) 


B = 


= ksn, 




(A3) 


C = 


= {akfi + bkc + ckB)n - 


= ken, 


(A4) 



where kfi, ka, fcs are integral quantities that can be com- 
puted from known spectra or directly measured (here 
they are obtained from experimental calibrations in pure 
gases at known temperatures). 

Relation (|A4[) also holds for a two-gas mixture and, 
if the emission of the two gases is decoupled, the total 
intensity will be the sum of the individual ones. Labeling 
the ambient gas as amb and the jet gas as jet, emission 
C can be written as the sum of two expressions of the 
kind [KM: 



(-^ — '-^amb + ^jet — kambnamb + kj^t^jet (A5) 
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and again using equation (|A4|) wc obtain 

ken = kambnamb + fcjet"jet- (A6) 

Dividing by the total numerical density n and introducing 
the concentrations Zamb ~ namb/n and Zjet = njet/n, 
equation (|A6|) becomes 



^C — ^amb^amb i f^jet^jet^ 



where 



•^jet 



akuarnb + bkaamb + cksa 

akiijet + bkajet + cksjet- 



(A7) 



(A8) 
(A9) 



Noting that the ratio between a pair of color intensities 
takes different values for different gases, it follows that 
two linear superpositions of the kind 



Ci =aiR + hiG + ciB 
C2 = 02^? + h2G + C2B 



(AlO) 
(All) 



will usually give a ratio C1/C2 that varies with the 
species concentration. For a given pair of gases, a suit- 
able choice of the six coefficients ai,bi,Ci {i = 1,2) will 
give a C1/C2 ratio that varies over the largest possible 
interval. This in turn leads to the concentration deter- 
mination, if Ci /C2 is rewritten using equations (jA4p and 
(!X7l) : 



L"! f^Cl f^ambi^amb ~r KjetiZjet 

L'2 "^02 kamb2 ^amb + "-^6*2 ^jet 



(A12) 



The concentrations of the two gases arc hnked by Zamb 
1 — Zjet'j writing C1/C2 = r, equation (|A12p becomes 






(A13) 



which can be solved with respect to Zjet, obtaining 



Zjet 



ib.'T 



^arnbi 



^jeti 



(k 



jet2 



f£n 



>.)r- 



(A14) 



This gives the concentration in the jet gas, for each pixel 
of the image, as a function of the ratio r of the two super- 
positions Ci , C2 read on the same pixel, provided that the 
four coefficients kambi,kamb2,kjeti,kjet2, defined through 
(jA8P - (jA9p . are known through calibration. 

Once the concentrations arc known, the density can 
easily be determined, as the total number density n may 
be obtained from relations such as (|A4[) . e.g. one of the 
foUowings relations can be used: 



Ci 


= kcin, 


kc. 


— f^ambl^amb ~r rCjetl^jet (,-^tOJ 


C2 


== kci^, 


kc2 


= kamb2Zamb + kjet2Zjet- (A16) 



Finally, the density is finally given by the equation 

P = n{Zambmamb + ZjetTrijet) (A17) 

where mamb,iTijet are the molecular masses of the gases 
in the ambient and in the jet, respectively. 
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